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Abstract. An elasto- plasticity model with coupled hardening variables of strain 
type is presented. In the theoretical framework of generalized associativity, the 
formulation of this model is based on the introduction of two hardening variables 
with a coupled evolution. Even if the corresponding hardening rules are linear, the 
stress-strain hardening evolution is non-linear. The numerical implementation by a 
standard return mapping algorithm is discussed and some numerical simulations of 
cyclic behaviour in the univariate case are presented. 



1 Introduction 

Starting from the analysis of the dislocation phenomenon in metallic materi- 
als, Zarka and Casier pQ and Kabhou et al. |2J proposed an elasto-plasticity 
model ("four-parameter model") where, in addition to the usual kinematic 
hardening internal variable, a second strain like internal variable was intro- 
duced. It plays a role in a modified definition of the von Mises criterion and 
its evolution, defined by linear flow rules, is coupled with the one of the kine- 
matic hardening variable. The resulting elasto-plastic model depends only on 
four parameters. Its non-linear hardening behavior was studied in [5j and a 
parameter identification method using essentially a cyclic uniaxial test was 
presented in [5] . In this note, the thermodynamic formulation of the classi- 
cal elasto-plastic model with linear kinematic and isotropic hardening is first 
recalled. Then, by using the same theoretical framework, a generalization 
of the four-parameter model is suggested, relying on the introduction of an 
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additional isotropic hardening variable. Finally, a return mapping implemen- 
tation of the generalized model is presented and some numerical simulations 
are briefly discussed. 



2 Thermodynamic formulation of a plasticity model 
with linear kinematic/isotropic hardening 

Under the assumption of isothermal infinitesimal transformations and of 
isotropic material, the hydrostatic and the deviatoric responses can be treated 
separately (see, among others, [7]). Hence, the free energy density & can be 
split into its spherical part Wh and its deviatoric part ^f d . To obtain linear 
state equations, tf^ and tf'd are assumed quadratic. Moreover, experimental 
results for metals show that permanent strain is only due to deviatoric slip. 
Hence, an elastic spherical behaviour is assumed, leading to the following 
definition : 

^ = \{x+ 2 -^jtr{ef = \k tr{ef (1) 

where e is the (small) strain tensor, A and fx are the Lame constants and K is 
the bulk modulus. Under the same assumptions, the deviatoric potential ^ 
must depend only on deviatoric state variables. The plastic flow is associated 
to the plastic strain e p , while the kinematic/isotropic hardening behaviour is 
introduced by the tensorial internal variable a and by the scalar variable p : 

$d = &d (e d ,e p ,a,p) = -t(e d - e p ) : (s d - e p ) + — a : a+—p 2 (2) 

where tr (a) = tr (e' p ) = and B, H > . The evolution of p will be related 
to the norm of e p . 

The state equation concerning the deviatoric stress tensor is easily derived: 
a d = ^ = 2^(e d -e p ) (3) 

0£ d 

and the thermodynamic forces associated to e p , a and p are defined by : 
cj d = -^ = 2^(e d -e p ) 

*=lt=Ba (4) 

One can notice that tr (e) — 0. The linearity of the hardening rules ([31)2-3 
follows from the quadratic form assumed for the last two terms in The 
second principle of thermodynamics can be written as follows [7] : 

(id ■ e d - W d > (5) 
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By using ((4|) in J5]), the Clausius Duhem inequality is obtained : 

a d : i p - X : a - R p > (6) 

In order to fulhl this inequality, a classical assumption is to impose that (s p , a, 
p) belongs to the subdifferential of a positive convex function cf> d , equal to zero 
in zero, called pseudo-potential. In such a case, the evolution of the internal 
variables is compatible with ([6]) [3]. 

The von Mises criterion corresponds to a special choice for the pseudo- 
potential (jA (a d ,~K., R) which is equal, in this case, to the indicator function 
I/<o of the elastic domain, or, to be more specific, of the set of (cr^X, i?) 
such that the so-called yielding function / is non-positive : 

/ = / (a d , X, R) = \\a d - Xj| - ^a y -R<0 (7) 

where ||-|| is the standard L2-norm. To impose that (i p ,a, p) belongs to the 
subdifferential of I/<o is equivalent to write : 

' i p = A a f = A 1Td ~ x 

d(Jd ||cr d -X|| 

<« = -A^=A^| r (8) 
[ p = -X§L = \ 

with the conditions A > , / < and A/ = . 

Equations © are called generalized associativity conditions or associative 
flow rules. The relations ([5]) yield in this case a = i p and A = p = \\i p \\. 
^From ©2 an d ©2 one obtains the Prager's linear kinematic hardening rule 
and a linear isotropic hardening rule : 

± = Be p , R = Hp (9) 



The coefficient A is strictly positive only if / = 0. In this case, its value can 
be derived from the so-called consistency condition / = 0, i.e. 

df . df • df • 
d^ d :(7d+ dX :X+ dR R = 

The introduction into the previous equation of the state equation ([3]) , as well 
as the thermodynamic force definitions ((4j)2— 3 and the normality ((8]), yield : 

14— 11=» 



Moreover, in a strain driven approach, Eq. (|I0|) has to be rewritten still using 
the state equation ([3|). As a result, by collecting A, one obtains 

n[J> 2u ^-^ + B^-^ + H^^ l + M+K Ikrf-XII " 
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where 7i (/) is zero when / < and equal to 1 for / = 0. The symbol (.) 
represents the MacCauley brackets. 



3 A generalization of the four-parameter model 

The linear hardening model discussed previously is used here to suggest a 
generalization of the 4-parameter model cited in the introduction. The tensor 
a into Eq. d2|) is replaced by a couple of tensors (cci,^)- As a result, the 
scalar constant B becomes a 2 x 2 symmetric positive definite matrix, denoted 
by B — [bij]. For sake of simplicity, only the thermodynamic potential is 
considered here and it is defined as: 

& d (e d ,E p ,ai,a 2 ,p) = y (e d - e p ) : (e d - e p ) + ^a T B a+|p 2 

where [i and H have the same meaning as before and a is the column vector 
defined as a = [a\; a^\. The state equation becomes : 

<T i = 1 r± = 2n(e d -f?) (11) 
oe d 

and the thermodynamic forces have the following form : 

Xl = ^ = 6 11 ai+&i 2 a2 (a d = 2^{e d -e p ) 

Y U\ h 21 OT { x = B a ( 12 ) 

X 2 = ssf = b 2i r «i + 622 «2 \ R-H v 

where X = [Xi; X2]. The Clausius -Duhem inequality becomes in this case : 
cr rf : e p - Xi : a x - X 2 : a 2 - R p > or a d : e p - X T a - R p > 



Moreover, the loading function / is defined as follows : 



/ = /(c r d,Xi,X 2 ,i?) = 




-R<0 (13) 



with p a positive scalar. One can remark that for p = and H = the 
standard von Mises criterion is derived while for p = 1 and H = the 
4-parameter model is retrieved. The flow rules are defined by a normality 
condition : 



(e p , -at, -d 2 , -p) £ d4>* d = <9I/< 



Therefore, the proposed model belongs to the framework of generalized asso- 
ciative plasticity [3] . The loading function (fl~3"|) can be rewritten as 
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/ = 3 (Y 1 ,Y 2 ) 



d<Jd Vll CT d- X lH 2 +P 2 H X 2l 

°i - ~ A ax7 _ A 



R < with Yi = a d - Xi and Y 2 



-X 2 . Hence 



x /|| - d -X 1 || 2 +p2||X 2 || 2 



< 8/ ; p 2 x, 



Vlkd-Xif+^HXa 



d = AVg 
P = A 



(14) 



It can be seen from (IT41 that &i = e p and p = A = ||d| = y ||ai|| + ||a 2 | 
Moreover, from (fl"4|) i_9 and (fT2|) 9_A one obtains the kinematic and isotropic 
hardening rules : 



Xi = bni p + 612 d 2 , X 2 = b 2 ii p + 6 22 d 2 , 
In [4] it was proved that B can be written as : 



R = H p 



B 



(A 



- rb 



r 2 6) — rb 



where the scalars A x and b are strictly positive and have the dimension 
of stresses. For r = 0, there is no coupling and B is diagonal, so that the 
dimcnsionless scalar r can be seen as a coupling factor in the evolutions of 
Xi and X 2 : 



Xi = e p + rb (ri p -a 2 ), X 2 = b(ri p - d 2 ), 



"X 2 = A r _ 



In the first two flow rules a recalling term appears, as in the non-linear 
kinematic hardening model of Frederich and Armstrong [5]. As before, the 
plastic multiplier can be explicitly computed by the consistency condition : 



A = «(/)■ 



Vg.B.Vg+H 

111 



> 0. 



4 Implementation and some numerical results 

In this section, a numerical implementation of the model is proposed. A 
standard return mapping algorithm is considered (see :6j). The formulation 
is explicitly described in the univariate case, but the tensorial generalization is 
straightforward. Let At n be the amplitude of the time step defined by t n and 
t n+1 and let d„ = [ai l7l , a 2j „,p„] T and X„ = [Xi >n , AT 2i „, R n ] T be the vectors 
collecting the internal variables and the corresponding thermodynamic forces. 
Moreover, let 
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be the global hardening modulus matrix. In a strain driven approach, knowing 
the value of all the variables at the time t n and the strain increment Ae n 
occurring during the time step t n — > t n +i, the numerical scheme computes 
the variables value at t n +i: 



(e n ,sl,a n ,a n ,X n J n ) + As n ^ ( £n+l i £ n +i j OLn+lt °n+l) ^n+lj /n+1 



The flow equations (Tl4ll define a first order differential system, which can 
be solved by the implicit Euler method. Therefore, the discrete form of the 
model evolution rules is (the notation 9 W / is equivalent to df/dw) : 

Rn+l)<0 



fn+l '■= \/(o"n+l— Xi, n+ i) 2 + p 2 (X 2 ,n + l) 2 — (p y 



'n+1 



'n+1 



-- E (e n+1 - ef l+1 ) ; 

e n = Aln+l da. fn+l 



X, 



„_ ] = D5„+i with E = M ^x+j^ 



a n +i — a n = — Ay n+ i 8%f n+ i 



A"f n+ i > 0, f n+ i < 0, Aj n+ i /, 



n+1 



0. 



An elastic predictor-plastic corrector algorithm is used to take into ac- 
count the Kuhn- Tucker conditions [6] (cf the last row). At every time step, 
in the first predictor phase it holds f n <0, an clastic behaviour is assumed 
and a trial value of f n +i , i.e. fj^h, is computed. If f^li < , then an elastic 
behaviour occurs, A"f n+ i has to be zero and no corrector phase is required. 
On the other hand, if f^li > , then plastic strains occur, the elastic pre- 
diction has to be corrected and Z\7„+i > has to be computed. This is done 
by a suitable return mapping algorithm, described below : 



i) Initialization 



ft, — u, t n+1 — fc n> u n+l — u„, /, 
ii) Check yield condition and evaluate residuals 



x(o) 



„(°) - 







(fc) 



:= £ 



R 



CO 

n+1 



+i " 
-p(fe) 



;(fe) 



-n+1 



— D 5 (fc) 

1 •— ^ a n+l 



' Jn+1 



if: f { n % < toh 



& 



+ ^ W 

t 7n+l 



R 



(fc) 

n+1 



daf 

< to^ 2 



f fn- (fe) ic (fe) ^ 

J {V-n+V A n+1 J 



then: EXIT 



iii) Elastic moduli and consistent tangent moduli 



(fc) 

n+1 



E 



D 



(fc) 



D 



(AS,)"' 



'n+1 

{Cnll+ln+ldl a .fn+l) 



7n+l<9^/„+l 



"n+l^/n+l ( D n+1 + 7n+l f n +l 



XX' 



iv) Increment of the consistency parameter 



^/n+l — 



Jn + l [OgJn+1 d y^Jn+l\ ^n+l^n+1 



[ftr/ffl ajt/Sl]'^!^?! ^x/n+lf 

v) Increments of plastic strain and internal variables 



CO 
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D 



71+1 



(fe) 



Ak) 

L Tl + l 



R 



(fe) 

n+l 



(fe) 
+ 1 



a f (fe) 

"oVn+l 
"xJn+1 



vi) Update state variables and consistency parameter 



_ P (fc+i) _ P (fc) 



j>(fc) 



'n+l 



-(fe+i) _(fe) . A -(k) (fe+i) (fe) . A a 



(fe) 
+i 



This procedure to determine Ay n+ i requires the computation, at each it- 
eration, of the gradient and the Hessian matrix of / . Other algorithmic 
approaches by-pass the need of the Hessian of / , but they are not considered 
here. 
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This implementation is used to obtain hysteresis loops in some particular 
cases. The values of the four parameters E, <j yi Aoo and b are the same as those 
used in [5] and correspond to the identified values of an Inconel alloy (E = 
205580 Mpa, a y = 1708,9 Mpa, = 35500 Mpa, b = 380700 Mpa). The 
value of the new parameter p is p = 1 and the values of r and H are indicated 
in the caption of each figure, /newline Fig. Q] frustrates the hysteresis loops 
obtained with an increasing amplitude strain history. The effect of the newly 
introduced isotropic hardening term is highlighted. Fig. [2] refers to a stress 
input history, with constant amplitude and non-zero mean. The plastic strain 
accumulation (ratchetting) and the shakedown phenomenon are modelled by 
changing only one parameter. The hysteresis loops are qualitatively similar 
to the ones of the non-linear kinematic hardening model of Armstrong and 
Frederick [8]. 



a) 
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Fig. 1. Hysteresis loops for an imposed history with increasing strain amplitude, 
a) r = 0.608, H = MPa b) r = 0.608, H = 6500 MPa. 



5 Conclusions 

A model with coupled hardening variables of strain type has been presented. 
It permits to take into account isotropic hardening and to have an elastic 
unloading path of varying length depending on the history of the loading. The 
simplicity of this model, which depends only on six parameters, seems to be 
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very attractive for structural modelling applications with ratchetting effects. 
To this aim, the proposed return mapping algorithm is a useful numerical 
tool, which allows numerical simulations to be performed in an effective way. 




2 4 6 8 
plastic strain (epsp) 

Fig. 2. Hysteresis loops for an imposed history with constant stress amplitude and 
non-zero mean stress, a) r = 0.608, H = MPa; b) r = 0.9, H = MPa. 
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